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ABSTRACT 



One method of reducing truncation error in numerical 
prediction is the use of more accurate finite-difference 
approximations. A relatively simple barotropic model, 
similar to the one currently in use by the Fleet Numerical 
Weather Facility (FNWF), Monterey, California, is employed 
as a basic program in testing several higher order finite- 
differencing schemes. Forecasts were computed to 24 and 
4$ hours with modifications of the basic model, and 
comparisons were made with the analyses at verifying time. 
It is the intent of this study to determine the effects of 
some higher order finite-difference approximations in this 
numerical prediction model. 

The author expresses his appreciation to Professor 
George J. Haltiner for his suggestions, guidance, and 
patience. Appreciation is also expressed to the personnel 
of FNWF for providing data and programming assistance. 
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I. Introduction. 

In any numerical method of weather prediction, 
regardless of the simplicity or complexity of the forecast 
model, one source of error which cannot be entirely 
eliminated is the so-called truncation error. This error 
is mathematical in nature and stems from the use of finite- 
difference approximations to the derivatives of functions 
such as those encountered in meteorology. Thompson [7] 
discusses the problem of truncation error and outlines 
two methods of reducing the magnitude of the error. The 
first method involves decreasing the grid interval thereby 
increasing the number of grid points for a given area. 

This procedure is not considered in this study. The second 
method utilizes the "quasi-analytic" property of the 
functions dealt with in numerical prediction which implies 
that the functions and their derivatives vary smoothly. 

If in reality these functions are analytic, uniformly 
convergent Taylor series can be written for any such 
functions. The inclusion of higher order terms of the 
Taylor expansion results in more accurate finite-difference 
approximations . 

It is the purpose of this study to determine the 
effects of some higher order finite-difference approximations 
in numerical weather prediction with the intention of 
obtaining a more accurate prognostication. 
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2 . 



The Basic Model 



The basic model used in this study is a simple model 
described by Haltiner [4] and similar to the 500-mb 
barotropic-type forecast currently in operational use by 
the Fleet Numerical Weather Facility (FNWF), Monterey, 
California [2] . This model was selected because of its 
relative simplicity, availability of data in a convenient 
form, and suitability in programming . It was supposed that 
a simple model would lend itself best to this type of study. 
In this model, three empirical terms are incorporated; 

(a) a "Helmholtz" term, (b) a "latitudinal" term, and 
(c) a "temperature" term. Although these terms are 
reported to be relatively small, they are included because 
they have been found to improve the results. 

The equation programmed for the basic model is of the 

form : 



Equation (1) is solved numerically by a sequential 
relaxation method. The input parameters are the 500-mb 
analysis and the thickness of the layer between 1000 mb and 
500 mb, with only the values north of 2N entering into the 
computations. A forward time step is used initially 
followed by centered-differences at one-hour intervals. At 
six-hour intervals, slight smoothing is applied and the 
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ellipticity criterion imposed, after which a forward time 
step is again initiated and the process repeated. 

The finite-difference approximations programmed in 
the basic forecast scheme can be found in most texts on 
numerical methods. Collatz [l] provides a comprehensive 
list of these formulas. Formulas used in the basic program 
are the 5-point approximation to the Laplacian 



V*Z„ - (z, + Z 2 *Z^Z,-4Z 0 ) 

and the 5-point approximation to the Jacobian 
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where the numerical subscripts refer to the grid position 

with the 5-point grid defined as 
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The time stepping formulas consist of the forward-difference 
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and the centered-difference 
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where the subscripts refer to the time at which the 
quantities are computed. 

To reduce the truncation error in the basic forecast 
program many higher order finite-difference formulas are 
available. Only a few of these formulas are considered 
in this investigation. To identify the modifications to 
the basic program, letter designations were given to the 
changes which are described in the following section. 

3. Modification B. 

The first modification to the basic program consists 
of substituting for equation (2) a 9-point approximation 
to the Laplacian 

V 2 Z. (6a , 

Equation (6a) can also be written symbolically as suggested 
by Forsythe and Wasow [3]: 

-/ 

/(, 

-I /6 -6o !(o -I H. 

IU 

-/ J (6b) 

The remainder, representing the residual between the 
finite-difference approximation and the actual value, 
consists of sixth-order derivatives and higher. 

4. Modification C. 

The same grid as presented in section 3 can be used to 
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